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I . INTRODUCTION 


If  one  is  interested  in  the  distribution  of  the  sum  of  two  random 
variables,  an  immediate  approach  would  be  to  use  Fourier  transforms. 

It  is  well  known  that  the  Fourier  transform  of  the  density  function  of 
the  sum  of  two  random  variables  is  equal  to  the  product  of  the  Fourier 
transforms  of  the  density  functions  of  the  respective  random  variables. 

To  determine  the  distribution  of  the  product  of  two  random  variables, 
an  analogous  approach,  namely,  the  application  of  Mellin  transforms,  will 
be  given  in  this  paper.  In  addition,  a recursive  relationship  will  be 
developed  to  evaluate  the  resulting  distribution. 


II.  DEVELOPMENT  OF  DISTRIBUTION 


(5) 


The  distribution  (8)  is  parameterized  by  an  integer  n.  We  will  now 
proceed  to  develop  a recursive  relationship  in  terms  of  the  parameter 
that  will  facilitate  the 

Define 

Pn(z)  = 

To  simplify  notation,  we 


JThis  problem  was  suggested  to  the  author  by  D.  Clark,  AMSAA. 

'y 

“Erdelyi,  A.,  et  al , Editor,  "Table  of  Integral  Transforms",  Bateman 
Manuscript  Project,  Vol . I,  McGraw-Hill  Book  Co.,  Inc.,  New  York,  1954 
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evaluation  of  the  distribution. 


r-,,3-2n  r2/z  . _ 


(n-l)!(n- 


(u)du. 


(9) 


will  consider 

• a 
'0 


I = f iTn  “k. (u)du. 
n Jn  1 


(10) 


From  the  identity  (see  [3]) 


K'(x)  = -K  (x)  + — K (x), 
J m+ 1 v ' x m v ' 


we  have,  where  from  (11)  (u)  = -Kg(u) , by  integration  by  parts 


I =-[*  u 

n Jo 


f u2n'2K' (u)du  = -a2n2K  (a)  + (2n-2)  f u2n_3K  (u)du. 
0 U J0  U 


From  (11),  where  m = -1  and  the  fact  that  K (x)  = K (x) , we  have 

mv  -nr  ' 

K0(x)  = -K’(x)  - x”1K1 (x) . 

If  we  substitute  (13)  into  (12)  and  integrate  by  parts,  we  obtain 
In  = -a2n'2KQ(a)  - (2n-2 )J  [u2n'3KJ(u)  + u2n'4Kj (u) ]du 


= -a2n_2K0(a)  - (2n-2)a2n’3K1(a)  + 


+ (2n-2) (2n-4) 


2n-4 

Jr 


(u)du.  (14) 


Therefore,  from  (9) 

Pn(z)  = ^ijrrn'-T)  ! [-(2^)2n~2K0(2/T)  - (2n-2)  (2/7) 2n~\  (2/?)  ♦ 


= P (z)  - 

n-lUJ  (n-1) 


+ (2n-2)  (2n-4)J  u2n'4KJ  (u)duj 
in^)!(n-2)l  V2^1  * <2"-2)  S <2/r>]  • 


for  n > 2,  where 


P2  = 1 - 2zK2(2/T). 


Abramowitz,  Milton  and  Stegun,  Irene  A.,  Editors,  "Handbook  of  Mathematical 
Functions  with  Formulas,  Graphs,  and  Mathematical  Tables",  National 
Bureau  of  Standards,  Applied  Mathematics  Series-55,  1964. 
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and  K^Cx) , m = 0,  1,  2,  are  mth  order  modified  Bessel  functions  of  the 
third  kind. 


IV.  CONCLUSION 


The  original  problem,  as  presented  to  the  author,  was  to  construct  a 
double-entry  table  of  the  distribution  over  the  variable  z and  the 
parameter  n.  A numerical  scheme  involved  an  approximation  of  a double 
integration  over  an  infinite  limit  had  been  programmed  for  BRLESC.  This 
procedure  was  estimated  to  take  an  average  of  five  minutes  of  computer 
time  to  evaluate  a single  point  in  the  double-entry  table.  Equations 
(15)  and  (16)  reduced  this  to  a table  look-up  and  a simple  calculation 
with  a hand  calculator. 
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